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ABSTRACT 

We report the detection of 3 additional Wolf-Rayet stars in the young cluster Westerlund 1. They were selected as emission-line star candidates 
based on 1 fim narrow-band imaging of the cluster carried out at OPD/LNA (Brazil), and then confirmed as Wolf-Rayet stars by K-band 
spectroscopy performed at the 4.1 m SOAR telescope (Chile). Together with previous works, this increases the population of Wolf-Rayet stars 
detected in the cluster to 22 members. Moreover, it is presented for the first time a K-band spectrum of the Luminous Blue Variable W243, 
which apparently implies in a higher temperature than that derived from optical spectra taken in 2003. The WC9 star WR-F was also observed, 
showing clear evidence of dust emission in the K-band. 
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1. Introduction 

Massive stars have an outstanding impact in the chemical and 
dynamical evolution of a galaxy, and certai nly played a key 
role i n the reionization of the early universe ( Loeb & Barkana 
2001). Therefore, the precise knowledge of how massive stars 
evolve in the HR diagram, and what are their physical param- 
eters during their lifetime, is crucial to understanding the star- 
formation history of the universe. Massive stars, their strong 
winds, their large flux of ionizing photons, and their spectac- 
ular death as supernovae, are also the main characters acting 
in harsh environments such as starburst galaxies, Wolf-Rayet 
(WR) galaxies and super star clusters. 

In the last decades, the starburst phenomenon and its 
associated massive stellar population has been analyzed 
in Local Group galaxies as the LMC dMelnickl 1 19851 



Walborn & Parker 1992; Brandl et al. 1996) and beyond ( Conti 
1991 ; Watson etalJll996l IWhitmore et al] Il999t iMelo et alJ 



2005). With the tremendous development of the infrared as- 
tronomy in the last decade, the efforts to observing massive 
stella r clust ers turned towards our Galaxy IFiger et al.l ll999. 
I2002U2005I) . where the recognition of a massive young cluster 
has been challenging due to the high reddening, uncertainties 
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in the distance, and large spatial areas needed to be covered by 
surveys. 



This is probably the reason why the young stellar clus- 
ter Westerlund 1 (heenceforth Wd 1) remained in the shad- 
ows for al most 15 yea r s, sin ce the pioneering spectroscopic 
survey of IWesterlun5 Jl987l) suggested the presence of a 
large number of early-type stars. However, after the discov- 
ery of a significant popul ation of WR stars in the cluster 
dClark & NegueraeldEjol . Wd 1 has been the target of in- 
tense observations. It was revealed that Wd 1 is very massive 
(~ 10 5 Mo), harboring a h uge population of evolved mas- 
sive stars ( Clark et al. 2005). The la st census counts for 19 
WR stars ( Negueruela & Clark 2005), a luminous blue vari- 
able (L BV) star dC lark & Negueruela 2004), 6 yellow hyper- 
giants ( Clark et al. 2005}), and a recently-discovered neutron 
star {Muno et a ll2006h . 



The goal of this work was to detect obscured emission-line 
stars in Wd 1 . The 1 yum narrow-band imaging of the cluster, 
which was used to select the candidates, and the follow-up K- 
band spectroscopy are described in Sect. [2 In Sect. [2]it is pre- 
sented the spectra of the newly-discovered WRs, line identifi- 
cations, and spectroscopic classification. It is shown in Sect.0] 
for the first time, the K-band spectra of the LBV W243 and of 
the WC9 star WR-F. The conclusions are presented in Sect. [3] 
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Table 1. Properties of the 1 //m narrow-band filters used in this 
work 



Filter 


Central Wavelength (jum) 


FWHM Oum) 


Contl 


0.9900 


0.0194 


Hen 


1.0137 


0.0088 


Cont2 


1.0655 


0.0192 


He i 


1.0825 


0.0093 



2. Observations 

2.1. Imaging 

We carried out near-infrared imaging of the cluster at the 
0.6 m telescope of the OPD/LNA (Brazil), during the nights 
of 2004 June 27-29. The detector used was a HAWAII 
HgCdTe 1024x1024 infrared array (CamlV), with a pixel size 
of 18.5//m and plate scale of 0.5" pixel -1 . The images were 
obtained with narrow-band filters centred approximately in 
the emission lines of Hen 1.0124/iin and Hei 1.0830/mi. The 
continuum was sampled with narrow-band filters centered in 
0.9900 /mi and 1.0655 pun. The photometric system and this 
technique to dete ct emission-lin e objects was partially de- 
scribed by Damin eli et alJ 41997). The other narrow-band fil- 
ters used in this work (Hei 1.0830/im and Cont2 at 1.0655 /mi) 
follow a similar procedure as described by those authors. In 
Tabled we summarize the properties of the set of filters used 
in this work. 

The observing strategy was done following the standard 
techniques of near-infrared imaging. For each filter, a mosaic 
of 4 images of 60s was obtained, dithered by 10" in a square 
pattern. The frames were median-combined using a rejection 
criterion to discard the stellar contribution, obtaining a mean 
sky image. This sky image was subtracted from each original 
image, which were subsequently divided by a normalized flat- 
field image. The instrumental magnitudes of th e stars were ex - 
tracted using the IRAF 1 package DAOPHOT JStetsonlll987ll . 
The analysis was restricted to objects with errors in the instru- 
mental magnitude less than 0.05 mag. 

Once the instrumental magnitudes of the stars were ob- 
tained, line-continuum versus continuum diagrams were con- 
structed. An absolute calibration of the magnitudes was not 
attempted, as the aim was only to obtain a difference in 
magnitude between line-continuum. In Figure ^ it is shown 
this diagram for the filters centred at Hen 1.0 124 /mi and 
Hei 1.0830/mi. With this procedure, one can separate stars 
without emission lines from stars which actually have them. 

It can clearly be seen in Fig. [J th at 14 of the 1 9 WRs pre- 
viously detected in Wd 1 b y Clark & Ne gueruelal {2002) and 
iNegueruela & Clarkl ^2005) could be isolated in at least one of 
the filters. Those stars were labeled following the designation 
of these authors, and, for clarity, hereafter it is used their identi- 
fication letter. Among the 5 known WRs that were not detected 
using the aforementioned technique, WR-J, WR-K, and WR-R 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



are located in crowded regions, while WR-N lies outside the 
observed field. Interestingly, we could not detect any emission 
from the WNVL star WR- S. This is inconsistent with the op- 
tical spectra published by Negueruela & Clark (2005) which 
show detectable He I emission lines. Such a behaviour can be 
explained if this object is an LBV evolving towards the red 
side of the HR diagram. Therefore, we suggest that WR-S is an 
LBV candidate. 

It is also shown in Fig. ^ the locii of the LBV W243, the 
YHG W12, and the sgB[e] W9. The object W243 can barely 
be detected as a Hei 1.0830/im emission-line object, which 
is consistent with the equivalent width of this feature of just 
5 A, measured in spectra taken in 2004 July and 2005 July 
(Groh et al. 2006, in prep.). For instance, a typical WR has 
EW He 1 1.0830 ~ 300 A. The YHG W12 could not be resolved 
into its components a and b due to the poor image quality of 
our dataset. Moreover, the star WR-J is blended with W12, 
which probably explains why this YHG was detected as a 
He 1 1.0830/im emission-line source in Fig. [2 Alternatively, 
W12a itself could be an LBV going on an excursion to the blue 
side of the HR diagram, thus developing He i emission. To con- 
firm this scenario, spectroscopic and photometric observations 
of this object at high spatial resolution are required along the 
next years. 

2.2. K-band spectroscopy 

Interestingly, new emission-line candidates in Wd 1 could be 
identified from the narrow-band imaging. The spectroscopic 
follow-up of the candidates identified in Fig.[Qwas carried out 
using the SOAR 4.1m telescope atop Cerro Pachon, Chile. 

The near-infrared spectrograph OSIRIS 2 was used to gather 
long-slit spectra of the candidates at the K-band, centered at 
2.14/mi, with R=3000. This setup provided a wavelength cov- 
erage from 1.98 to 2.32 /im, which allow to access the most 
important diagnostic lines in the K-band. The spectra were 
obtained in 5 different positions of the detector, each one 
dithered by 5" from the previous, in order to remove bad pix- 
els and other detector features. The sky image was obtained 
by median-combining these dithered images with a suitable re- 
jection algorithm to remove the stellar contribution. The mean 
sky frame was subtracted from each original image, which was 
then divided by a normalized flat-field. The extraction of one- 
dimensional spectra was done using usual IRAF routines. The 
wavelength calibration was performed using the sky OH lines 
JOhva & Orig lialfl992). Telluric features were removed by di- 
viding the extracted spectra by the spectrum of a standard B 
star (HD 159402). The stellar Bry absorption feature present 
in HD 159402 was removed previously by interpolating the red 
and blue adjacent continuum. The individual spectra taken at 
different positions in the array were then combined to a final 

2 OSIRIS (Ohio State Infrared Imager and Spectrograph) is a col- 
laborative project between Ohio State University and CTIO. Osiris 
was developed through NSF grants AST 90-16112 and AST 92- 
18449. OSIRIS is described in the instrument manuals found on the 
CTIO Web site at http://www.ctio.noao.edu See also iDepov et alJ 
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Fig. 1. Line-continuum versus continuum diagram for the 
narrow-band filters used in this work. The units of both 
axis are arbitrary magnitudes. Upper panel: filter cen- 
tered at Hen 1.0 124 yum. Lower panel: filter centered at 
He 1 1.0830jL/m. It is noticeable, in both filters, the presence of 
new emission-line candidat es other than the WRs presented by 
iNeeueruela & ClarkllEo OS): the latter are identified by their re- 
spective identification letters. 

spectrum. For each object, the final spectrum was continuum- 
normalized with a low-order Legendre polynomial function. 

3. Newly-identified Wolf-Rayet stars 

In Figure [2] it is displayed the K-band spectra of the newly- 
discovered Wolf-Rayet stars in Wd 1 and the identification of 
the strongest emission lines. The stars were class ified based 
on the K-band spectral atlas of WRs published bv iFiger et alJ 
lll997l) . It can be readily seen that all of them are WN stars. In 
Table[2]it is presented the basic properties of the WRs observed 
with the SOAR telescope, which were named adding the prefix 
GDTB. A K-band finding chart is provided in Figure[3] 

The spectrum of GDTB # 1 has a lower S/N than the 
other WRs shown here due to its serendipitous detection 
while observing another candidate. The presence of strong 
Nv 2.0985/mi is ruled out, which implies in a subtype later 
than WN4. Strong lines of He i are also absent, which yields 
a subtype earlier than WN8. Hence, GDTB # 1 is tentatively 
classified as a WN5-7. 

The object GDTB #2 is the earliest WN of the sample, 
and probably one of the earliest WNs det ected in Wd 1, ri- 
valing with WR-A and WR-Q reported by Negueruela & Clark 
(2005) (see also Fig. The relatively strong emission of 
N v 2.0985 pm reveals its early type nature, together with the 
weakness of the blend at 2. 166 pun (H i+He i+He n) compared 
to Hen 2.1885/vm. The precise determination of the physical 
parameters of GDTB # 2, as of the other WR stars present in 
Wd 1, relies on a quantitative spectroscopic analysis, which is 
beyond the scope of this paper. 

The object GDTB #3 shows a typical WN7 spectrum. 
While the presence of He 1 2.058 1 pm is doubtful, H i is clearly 



Table 2. Basic data of the newly-discovered Wolf-Rayet stars 
in Westerlund 1. Coordinates and Ks magnitudes were ob- 
tained from images taken at NTT/ESO (M. Teodoro et al., in 
prep.). 



Object 


a (J2000) 


6 (J2000) 


K s 


Spectral type 


GDTB #1 


16h47m06.6s 


-45°50'38".6 


9.19 


WN5-7 


GDTB #2 


16h47ml4.2s 


-45°48'31".4 


9.99 


WN4-5 


GDTB #3 


16h47m07.6s 


-45°49'21".7 


9.70 


WN7 




wavelength l(lml 



Fig. 2. Continuum-normalized spectra of the 3 new WRs de- 
tected in this work. For the sake of clarity, the star GDTB # 1 
was shifted up by 2 dex and GDTB #2 by 1 dex. It can be 
noticed that they are all of the WN subtype. 

present in the blending at 2. 166 /mi. The other emission lines 
present in the spectrum are mainly due to He n transitions. Thi s 
object was detected as an x-ray source bv lSkinner et alJ J2OO6). 

4. Other emission-line objects 

As a part of the spectroscopic follow-up of the WR candidates, 
it was also obtained K-band spectra of other peculiar emission- 
line stars in Wd 1. They are shown in Figure |4] 

It is presented in Fig. 0] the K-band spectrum of the WC9 
star WR-F, which has an optical spectru m dominated by strong 
C 11 and C m lines JClark & Negueruelal2002l) . However, in its 
near-infrared spectrum, WR-F shows very weak He 1, He n, C 111 
and C iv features. This is consistent with the presence of warm 
dust (T ~ 10 00 A") emitting on the K-band, which veils the 
emission lines llFigeret alJll997l) . Moreover, the recent claims 
for t he presence of a puta tive companion star in dusty WC9 
stars (Williams et al. 2005) would yield another source of con- 
tinuum emission, which is compatible with the low-strength na- 
ture of the emission lines of WR-F. It was also recently detected 
as a moderately bright x-ray source in a Chandra observation 
of Wd 1 and a comparison of its X-ray properties with other 
WC stars suggests th at WR-F is likely a colliding-wind binary 
JSkinneretalJ200.5lE006li. 

A K-band spectrum of the luminous A-supergiant 
W243, classified as an LBV with spectral type A2I by 
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Fig. 3. K-band finding chart for the new WRs discovered in 
Wd 1. This image was obtained with NTT/ESO. North is up 
and east is to the left. The field covers a region of 4.5'x4.5'. 
The inset shows a box of 15"xl5" around the star W16. 

IClark & N egueruela (2004), was also gathered at the SOAR 
telescope. It shows strong, narrow emission lines (unresolved at 
our resolution of 100 km s~') of He 1 2.0581 /im and Bry. Weak 
emissions of Nai AA 2.2056-2.2084yum and MgnAA 2.1369- 
2.1432yum doublets can also be seen. The spectral morphol- 
ogy at the K-band i s very similar to other LB Vs, such as AG 
Car, LBV 1806-2 feikenberrv et alJl2004 and the Pistol Star 
jFigeret alJl998h . 

In Fig. |4] it is also shown a comparison between the K- 
band spectrum o f W243 and the YHG IRC 10+420 presented 
bv lHanson et alJ £l996). It can be noticed that IRC 10+420 
has a stronger emission of Nai than W243. However, W243 
has strong He 1 2.058 1 /im and Bry emission, which are both 
very weak in IRC 10+420. Hence, the K-band spectrum of 
W243 suggests a higher temperature than that of a typical YHG 
such as IRC 10 +420 (J, ff <10kK) or tha n that derived from 
optical spectra Jdark & Negueruelall2004l) . Indeed, this sug- 
gestion is supported by the comparison of the K-band spec- 
trum of W243 with supergiants that have temperatures ob- 
tained via non-LTE mode ls, such as AFGL 2298 (T e ff =12- 
15kK. IClarketal ll2003hh and G26.47+0.02 (T e// =17kK, 
IClark et alJl2003al) . Long-term spectroscopic and photometric 
monitoring of W243, especially in the near-infrared, is crucial 
to determine its current evolutionary phase and to estimate its 
physical parameters. 

5. Concluding remarks 

In this work it was reported the discovery of 3 additional WR 
stars in the young starburst cluster Westerlund 1 . Their K-band 
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Fig. 4. Continuum-normalized spectra of the other emission- 
line stars observed in Wd 1 . From top to bottom, it is shown 
the K-band spectrum of the WC9 star WR-F (shifted up by 0.6 
dex), of W243 (LBV) and a comparison between the spectra of 
IRC 10+420 (kindly provided by Dr. M. Hanson) and W 243 
(degraded to match the resolution of IRC 10+420). The formers 
were shifted down by 0.3 dex. 

spectra show that they belong to the WN subtype, in the range 
WN4-7. This i ncreases the po pulation of WRs detected within 
Wd 1 from 19 jNegueruela & ClarkbOoH) to 22 members. This 
work suggests that there is still a number of WRs to be discov- 
ered in this obscured super stellar cluster, since the narrow- 
band imaging was performed with a modest-sized telescope 
(0.6 m), reaching a limiting magnitude of K st 10.5 mag. 

It was also presented, for the first time, a K-band spectrum 
of a WC9 star in Wd 1 (WR-F), which clearly show evidence 
of continuum emission due to warm dust. This reduces the 
strength of the emission-lines in the observed spectrum, and 
potentially makes it difficult to detect another putative dusty- 
makers WC stars in the cluster through emission-line surveys 
in the K-band. Hence, a number of late WC stars may still be 
waiting to be discovered in Wd 1 . 

The complete census of early-type stars in Wd 1 is highly 
desirable, as it provides an unique "lab" to study a coeval, 
chemically homogeneous, large population of massive stars in 
the Galaxy. Of particular interest is the precise determination of 
the ratio WC/WN and WR/O, since they are dependent of the 
metalicity of the environment ( Mevnet & Maeder 2005). Also, 
the current rate of binarity is highly desirable to be known, es- 
pecially among the massive members of the cluster. This could 
give insights about the preferred mode of massive star forma- 
tion, and about the frequency of stellar mergers in such mas- 
sive stellar clusters. This can likely be the case of the peculiar 
sgB[e] W9, as suggested by Clar k et alJj2005l) . 

Therefore, a multi-wavelength monitoring campaign of the 
massive members of Wd 1 is essential over the next years. 
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